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C3 1 ABSTRACT 



. Context. Flat-spectrum radio sources often show a core-jet structure on pc-scales. Individual jet components reveal predominantly 

outward directed motion. For the BL Lac object S5 0716+714 conflicting apparent velocities have been reported in the literature. This 
object is an intra-day variable source and suited to investigate a possible correlation between kinematic properties and flux-density 
variability on different timescales. 

Aims. We study the kinematics in the pc-scale jet of S5 0716+714 to determine the apparent speeds of the jet components based on 
a much improved data set. In addition, we search for correlations between the radio flux-density light curves and the morphological 
changes detected along the VLBI jet. 
Methods. We (re-)analyze 50 VLBI observations obtained with the VLBA at 5 different frequencies (5 - 43 GHz) between 1992.73 
and 2006.32. The data have been parameterized using circular Gaussian components. We analyze the jet component motion in detail 
taking care not only to account for motion in the radial but also in the orthogonal direction. We study the evolution of the jet ridge 
line and investigate the spectral properties of the individual components. We search for correlations between radio band light curves 
and the kinematic properties of the jet components. 
Results. We present an alternative kinematic scenario for jet component motion in S5 0716+714. We present evidence for the apparent 
stationarity of jet components (with regard to their core separation) with time. Jet components, however, do seem to move significantly 
non-radially with regard to their position angle and in a direction perpendicular to the major axis of the jet. We discuss a possible 
correlation between the long-term radio flux-density variability and apparent jet component motions. 

Conclusions. In S5 0716+714 an alternative motion scenario is proposed. With regard to the core separation, rather stationary 
components can fit the VLBI observations well. A new model to explain the observed motion with regard to the position angle 
is required. Based on the correlation between the longterm radio flux-density variability and the position angle evolution of a jet 
component, we conclude that a geometric contribution to the origin of the long-term variability might not be negligible. Subluminal 
motion has been reported for most of the TeV blazars. Our analysis also confirms this finding for the case of S5 0716+714. This result 
increases the number of TeV blazars showing apparent subluminal motion to 7. 
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1 . Introduction optical (Gupta et al., 2008) and X-ray light-curves (Gupta et al., 

■ ™, „ T T ™ • „r ^ , , , , -r „ xt-, 2009). Apparently correlated optical/radio quasi-periods of ~1 

X • ^ ^ ° bj f S5 , 07 16+7 14 S dShlft °; 3 1 ±0 -? 8 ' N ? S ° n 6t day lasting over several days have been observed by Quirrenbach 

b ' u u } u flux -. densl jy v f ablh ty at all wavelength ranges gt aJ {ml) Q ^ aJ (2m) a h iod of „ 

; that have been investigated so far (e g. Wagner & Witzel 1995). 3 q ± 3 and G ta et al (2009) find quasi .p er iodic os- 

It is still the only Active Galactic Nuclei (AGN) for which a si- cmtions (QPOs) of -30 minutes. As a possible origin of the 

multaneous change in the mode of the radio/optical variability x QPOs, shocks propagating along the jet have been dis- 

has been observed on intra-day time scales (Quirrenbach et al. cussed Mogt recentl S5 0716+714 has been detected in the 

1991) and is thus so fai- the best candidate for an intrinsic origin XeV ime b the MAGIC collaboration (Teshima et al. 2008, 

of intra-day variability (IDV, Wagner & Witzel 1995, and refer- Anderhub et al 2009) 

ences therein). S5 0716+714 has been investigated with Very Long Baseline 

Several authors claim the detection of periodicity or nearly pe- Interferometr (VLB I) at various radio frequencies with differ- 

riodic variations in radio and optical (Quirrenbach et al. 1991), ^ mays Qn parsec . scales it reveals a North-South directed jet 

* Change of name: Veronika Meyer changed her name to Veronika with several embedded components (e.g., Eckart et al. 1987). 

Kam by marriage On kpc-scales, the jet is more East- West oriented with some ev- 

** Member of the International Max Planck Research School (IMPRS) idence for halo-emission (e.g., Antonucci et al. 1986). A rather 

for Radio and Infrared Astronomy at the Universities of Bonn and large range in proper motions (0.05 mas/yr - 1.11 mas/yr) has 

Cologne been reported by different authors based on the investigation of 
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Table 1. Details of the observations presented in this paper (in epochs marked by a star, S5 0716+714 served as calibrator for 
observations targeted on NRAO 150). 



Epoch 


v [GHz] 


Array/Survey 


Re 


erence 


1992.73 


5.0 


VLBA CJF-Survey 


Britzen et al. 


(2007) 


1992.85 


22.2 


Global array 


Bach et al. 


(2005) 


1993.71 


22.2 


Global array 


Bach et al. 


(2005) 


1994.21 


8.4 


VLBA 


Bach et al. 


(2005) 


1994.67 


15.3 


VLBA MOJAVE/2cm Survey 


Kellermann et al. (1998); Zensus et al. 


(2002) 


1994.70 


5.0 


VLBA CJF-Survey 


Britzen et al. 


(2007) 


1995.15 


22.2 


VLBA 


Jorstad et al. 


(2001) 


1995.47 


22.2 


VLBA 


Jorstad et al. 


(2001) 


1995.65 


8.4, 22.2 


VLBA 


Bach et al. 


(2005) 


1996.34 


22.2 


VLBA 


Jorstad et al. 


(2001) 


1996.53 


15.3 


VLBA MOJAVE/2cm Survey 


Kellermann et al. (1998); Zensus et al. 


(2002) 


1996.63 


5.0 


VLBA/CJF-Survey 


Britzen et al. 


(2007) 


1996.82 


15.3 


VLBA MOJAVE/2cm Survey 


Kellermann et al. (1998); Zensus et al. 


(2002) 


1996.90 


22.2 


VLBA 


Jorstad et al. 


(2001) 


1997.03 


8.4 


VLBA 


Fey et al. 


(2000) 


1997.58 


22.2 


VLBA 


Jorstad et al. 


(2001) 


1997.93 


8.4 


VLBA 


Ros et al. 


(2001) 


1999.41 


8.4 


VLBA 


Ros et al. 


(2001) 


1999.55 


15.3 


VLBA MOJAVE/2cm Survey 


Kellermann et al. (1998); Zensus et al. 


(2002) 


1999.89 


5.0 


VLBA/CJF-Survey 


Britzen et al. 


(2007) 


2000.82 


5.0 


VLBA+Effelsberg 


Bach et al. 


(2005) 


2001.17 


15.3 


VLBA MOJAVE/2cm Survey 


Kellermann et al. (1998); Zensus et al. 


(2002) 


2002.48 


8.4, 22.2, 43.2 


VLBA* 


Agudo et al. 


(2007) 


2003.20 


8.4, 22.2, 43.2 


VLBA* 


Agudo et al. 


(2007) 


2003.49 


8.4, 22.2, 43.2 


VLBA* 


Agudo et al. 


(2007) 


2003.65 


15.3 


VLBA MOJAVE/2cm Survey 


Kellermann et al. 


(2004) 


2003.88 


22.2, 43.2 


VLBA 


Agudo et al. (A&A, 2008 ir 


l prep.) 


2004.44 


15.3 


VLBA MOJAVE/2cm Survey 


Kellermann et al. 


(2004) 


2004.60 


15.3,22.2, 43.2 


VLBA* 


Agudo et al. 


(2007) 


2004.79 


15.3 


VLBA MOJAVE/2cm Survey 


Kellermann et al. 


(2004) 


2004.97 


15.3,22.2, 43.2 


VLBA* 


Agudo et al. 


(2007) 


2005.01 


15.3 


VLBA MOJAVE/2cm Survey 


Kellermann et al. 


(2004) 


2005.09 


15.3,22.2, 43.2 


VLBA* 


Agudo et al. 


(2007) 


2005.70 


15.3 


VLBA MOJAVE/2cm Survey 


Kellermann et al. 


(2004) 


2005.72 


15.3 


VLBA MOJAVE/2cm Survey 


Kellermann et al. 


(2004) 


2006.32 


15.3 


VLBA MOJAVE/2cm Survey 


Kellermann et al. 


(2004) 



an increasing core separation with time modeled for jet com- 
ponents in this source (e.g., Eckart et al. 1987, Witzel et al. 
1988, Schalinski et al. 1992, Gabuzda et al. 1998, Jorstad et al. 
2001, Perez-Torres et al. 2004, Kellermann et al. 2004, Bach et 
al. 2005, Rastorgueva et al. 2008). Based on these results it is 
not clear, whether S5 07 16+7 14 reveals slow apparent velocities 
typical for BL Lac objects (e.g., Britzen et al. 2007) or faster 
apparent velocities more typical for quasars (e.g., Kellermann et 
al. 2004). 

Based on the now available improved database with multi- 
frequency observations at 5 different frequencies, covering a 
time span of ~14 years, we can address the question concern- 
ing the kinematics in the pc-scale jet of S5 0716+714 in more 
detail. This source belongs to the most rapidly variable radio 
sources (e.g., Wagner & Witzel 1995) and we can - in addition 
to the kinematics - investigate whether these sources show evi- 
dence for a possible correlation between morphological changes 
and flux-density variability. 

In an earlier publication of the results relying on the 15 GHz data 
only, we presented a description of preliminary results (Britzen 
et al. 2006). A detailed description of the data reduction, analy- 
sis, and investigation is presented in Meyer (2007). 
The current paper is organized as follows: we first give an 
overview of the data that have been (re-)investigated; we 



then describe the model-fitting results and discuss the model- 
component kinematics; we present and discuss a new motion 
scenario; and we investigate a possible correlation between the 
total flux-density data and the VLBI kinematics, and concentrate 
in the discussion on the implications and resulting consequences. 
Throughout this paper we use a Hubble constant of Ho = 71km 
s^'Mpc -1 . 

2. The observations, data reduction and analysis 

2.1. The observations 

S5 0716+714 has been observed in 34 VLBA observations and 
in 2 Global array VLBI observations between 1992.73 and 
2006.32. These observations have been performed at 5.0, 8.4, 
15.3, 22.2, and 43.2 GHz. 12 of the VLBA epochs have been ob- 
tained by the MOJAVE / 2 cm- VLBA survey group (MOJAVE, 
e.g., Lister & Homan 2005, 2cm Survey, e.g., Kellermann et al. 
2004). 4 observations were taken from the CJF-survey (Britzen 
et al. 2007). S5 0716+714 was used as calibrator source for 
observations targeted at NRAO 150 (Agudo et al. 2007) in 6 
epochs. Part of the data presented here have already been ana- 
lyzed and discussed by Bach et al. (2005). All of the above men- 
tioned data-sets have been re-analyzed by us. 
In Table [T]we summarize the observations presented in this pa- 
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per. In column (1) we list the epoch, in column (2) the frequency 
of the observation, column (3) lists the array and survey-type, 
and in column (4) we note references to the original publication 
of the individual data sets. We refer to the references in column 
(4) for more details concerning the observations and data reduc- 
tion. 
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Fig. 1. The jet ridge line (as denned by all component positions per epoch) is shown at 15 GHz in (a)-(d) and at 22 GHz in (e)-(h). 
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2.2. Model-fitting 

Model-fits have been carried out using the diftnap -package (e.g., 
Lovell 2000) and starting with a point-like model in all of the ob- 
servations. The final number of jet components necessary to fit 
the data adequately was reached when adding another jet compo- 
nent did not lead to a significant improvement in the value of chi- 
squared. We aimed at obtaining the optimum number of compo- 
nents for each individual data set. The quality of the observations 
can differ from epoch to epoch and, in addition, the flux-density 
of the source - in particular in this case of 0716+714 - is known 
to vary between epochs. A brighter core component can obscure 
part of the jet structure, and not all jet components might be vis- 
ible in this particular epoch. Vice versa, in the case of a faint 
core component when more jet components might be detectable. 
Thus, the optimum number of Gaussian components that fits the 
data best, differs between epochs. The uncertainties of the model 
component parameters have been determined by comparing the 
parameter ranges obtained by performing model-fits with dif- 
ferent numbers of model components (±1 component). Where 
this was not possible, we used for the core separation r and the 
position angle an uncertainty of 5 %. The uncertainties of the 
flux-density increase with increasing distance from the core. We 
accounted for this by adjusting the uncertainties accordingly: r< 
lmas: 5 %, 1< r 2mas: 10 - 20 %, and for r>2mas : 20%. We 
list all parameters of all re-analyzed epochs in Table [5] We list 
the frequency, epoch of observation, flux-density, core separa- 
tion, position angle, axis size, and component identification. In 
Figs. E] [8] 13 UU1 E] [HI and [13] we show the maps with the 
Gaussian components superimposed. 

2.3. Component identification 

The component identification was carried out based on the as- 
sumption that changes of the flux density, core separation, and 
position angle of the modeled jet components should be small on 
time scales between adjacent epochs. We identified the compo- 
nents across the epochs by taking all model-fit parameters into 
account. We show an example for the component identification 
in Fig. |2] The identified model components are marked by let- 
ters which are used throughout the paper for the identification 
scenario presented here. We are aware of the problems involved 
when tracing and identifying components through the epochs. A 
detailed discussion of the problems and uncertainties involved 
has been published by us in Britzen et al. (2008) with regard to 
the analysis of the kinematics of the CJF-survey (for 293 AGN). 
Whereas in most VLBI data analyses the small number of obser- 
vations hampers a proper component identification the situation 
is much better in the case of 0716+714 with a comparably large 
number of data sets. However, 07 16+7 14 is a rather faint AGN in 
the cm-regime, and the jet is even fainter. Thus the modeling of 
jet components is difficult, and the identification of components 
is even more difficult. Although a large number of VLBA obser- 
vations (50) is available for the analysis, the data are -as usual in 
VLB analysis - unevenly sampled. In addition to this fact, VLBI 
data are usually of non-uniform quality. The identification across 
the epochs has to take differences in the number of data points, 
uv-coverage, data quality and resolution of the data sets into ac- 
count. To prevent a potentially large systematic error from aris- 
ing due to the incorrect cross-identification of moving features 
from epoch to epoch, we chose to adopt the simplest scheme 
when identifying the jet-features in the jet of 0716+714. We used 
all available model-fit parameters (flux-density, core separation, 
position angle, size of axis) to find the best suited set of cross- 
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identifications. The result presented here is meant to provide the 
most robust and simplest identification based on the largest VLB 
data set currently available for this source. Stroboscopic effects 
can occurr in very fast sources where too few data are available. 
For this source we present here the analysis of the largest data set 
so far available and hope we can exclude stroboscopic effects. 



Clean I map. Array: BEFGHKLMNOPSTR 
0716+714 at 4.978 GHz 2000 Oct 04 




core 
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Right Ascension (mas) 
Map center: RA: 07 21 53.448, Dec: +71 20 36.373 (2000.0) 
Map peck: 0.466 Jy/beam 

Ccntours %: -0.05 0.05 0.1 0.2 0.4 0.8 1.6 3.2 6.4 

Contours X: 12.8 25.6 51.2 

Boom FWHM: 0.988 x 0.878 (mas) ot -37.3° 

Fig. 2. A map of S5 0716+714 with Gaussian model-fit compo- 
nents superimposed. Letters denote individual components. 



3. Results 

In the following sections, we present the results of our model- 
fitting analysis. S5 0716+714 reveals a core-jet structure with 
a rather faint jet. Tracing the jet components across the epochs 
allows us to investigate the kinematics in the pc-scale jet. In ad- 
dition to tracing the components, we investigate the jet ridge line 
as a whole. This does not require any component identification 
and is thus independent of any identification scenarios. We first 
describe this evolution of the jet ridge line, and then we concen- 
trate on the core separation behavior of the individual compo- 
nents as a function of time. In the next subsection we compare 
our identification scenario - based on the observed core separa- 
tion/time relation - with results presented by other authors. We 
then discuss the position angle changes, calculate the apparent 
velocities in the standard manner (core separation/time) and tak- 
ing all the observed motions into account (XY). We compare the 
derived values with the results presented by other authors. We 
then study the frequency dependence of the individual compo- 
nents and finally search for correlations between the kinematics 
and the flux-density evolution. 




Fig. 3. The core separation as a function of time is shown for all frequencies investigated: (a) at 5 GHz, (b) at 8 GHz, (c) at 15 GHz, 
(d) at 22 GHz, and (e) at 43 GHz. We show the new identification scenario described in detail in this paper (filled symbols, linear 
regression is shown as solid line). In addition we show - where possible - the previous jet component identifications adopted by 
Bach et al. (2005) (open symbols, linear regression as dotted line). 



3. 1. Evolution of the jet ridge line 

In Fig. Q] we show the evolution of the jet ridge line (defined 
by the positions in core separation and position angle for all 
the components per epoch) with time at 15 GHz (a-d) and at 
22 GHz (e-h). Obviously the components define rather different 
paths between adjacent epochs. The difference in position an- 



gle is larger than the uncertainties of the individual components. 
Based on these plots, it seems that the individual components do 
not follow a well-defined path, instead, the ridge line changes 
between adjacent epochs. We are currently investigating much 
denser sampled data (one epoch per day over one week of obser- 
vations) for S5 0716+714 to investigate this effect in more detail 
and to determine whether this evolution is periodic. 
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linear regression 


data 


v [GHz] 


time span 


2 


number of data points 


'stationary scenario' 


Meyer 2007 


22 


1995.15-1996.90 


0.015 


10 


Bach et al. (2005) 


Bach et al. (2005) 




1995.15-1996.90 


0.016 


29 



Table 4. Apparent velocities taken from Bach et al. (2005) 



Id. 


y[maslyr\ 


Aw [c] 


CI 


0.86±0.13 


16.13+2.36 


C2 


0.63±0.10 


11.87+1.82 


C3 


0.62+0.04 


11.59+0.76 


C4 


0.44+0.01 


8.27+0.24 


C5 


0.43+0.02 


8.16+0.29 


C6 


0.41+0.02 


7.80+0.33 


C7 


0.37+0.01 


6.89+0.19 


C8 


0.32+0.01 


6.04+0.20 


C9 


0.26+0.01 


4.98+0.27 


C10 


0.24+0.01 


4.52+0.45 


ell 


0.29+0.02 


5.48+0.55 



3.2. Core separation as function of time 

In Fig. Oa)-(e) we show the core separation as a function of 
time for all model-fitted components for all data sets from 5 

- 43 GHz. The components are shown with filled symbols. 
Superimposed are linear regressions (solid line) based on our 
component identification. In addition to our own results, the re- 
sults published by Bach et al. (2005) for a subset of the data 
presented here are shown as well. The component identification 
across the observing epochs by Bach et al. is based on a com- 
bined data set composed of all data at all frequencies (5, 8, 15, 
and 22 GHz). The data by Bach et al. are indicated by open sym- 
bols, the linear regression by Bach et al. is shown with dotted 
lines. 

The larger amount of data analyzed in this paper leads to a dif- 
ferent component identification compared to Bach et al. We find 
that components tend to remain at similar core separations over 
the entire observing period. This behavior is found at all five fre- 
quencies. 

3.2.1 . Comparison of goodness of fit 

In the previous section we introduced a new scenario of jet com- 
ponent identification for S5 0716+714. Based on this core sepa- 
ration/time relation we obtain different apparent velocities com- 
pared to Bach et al. These apparent velocities will be discussed 
in section 3.4. Here we compare how good the assumed iden- 
tification (outwards moving components in the case of Bach et 
al. (2005) and rather stationary components as assumed in this 
paper) fits the distribution of the observed data. We therefore 
compare the reduced chi-squared values obtained when fitting 
these two scenarios to the data sets with the most intense sam- 
pling in order to have the most reliable data sets for a meaningful 
analysis: 

- the "stationary scenario" for the 22 GHz data by Meyer et al. 
(2007) between 1995.15 and 1996.90 

- the "fast scenario" by Bach et al. (2005) for the data from Bach 
et al. (2005) between 1995.15 and 1996.90 

In Table [2] we list in the first column the kind of scenario we 
apply, the second column identifies the data reduction, the third 



column lists the frequency, the fourth column lists the time span 
covered by these observations, the fifth column gives the reduced 
chi-squared value and the last column contains the number of 
data points that entered this study. We do not list the frequency 
for the data of the Bach et al. paper, since data from all frequen- 
cies were taken together to identify individual components; thus, 
all frequencies enter the reduced chi-squared value test. This ex- 
plains why the number of data points is significantly higher com- 
pared to the data set from Meyer (2007). Based on the assumed 
scenario ("fast" or "stationary") we apply a linear regression to 
the data and measure the quality of the fit by means of the chi- 
squared values. The;f 2 values were calculated as follows 

2 ^(observed values - expected values) 2 
^ expected values 



number of data - 1 

xl eA is the reduced^- 2 value. 

The resulting values are listed in Table|2] We obtain a compa- 
rable (insignificantly better) chi-squared value for the "station- 
ary" scenario-based linear regression for the 22 GHz data com- 
pared to the Bach et al. "fast scenario" for the same time range. 
We conclude that our scenario of rather stationary components 
fits the observed data well, and the test using the time spans with 
the most intensive sampling demonstrates quantitatively that the 
fit is at least comparable to, if not better than, the "fast motion 
scenario". 

A general comparison of the two scenarios is complicated 
by several factors: 

- Bach et al. (2005) determined the linear regression for 
individual components based on the multi-frequency dataset. 
All component parameters obtained at all investigated fre- 
quencies have been taken together to determine individual 
component motions. 

We investigated each frequency individually. 

- The time span covered by Bach et al. is significantly shorter 
compared to the time span that the data presented by us 
cover. 

- The number of components in Bach et al. is much larger 
compared to the number of components which our identifi- 
cation requires. 



3.3. Position angle as function of time 

All jet components show significant changes with regard to the 
position angle. We show the motion of the inner four jet compo- 
nents at 15 GHz in rectangular coordinates in Fig. |4] We used a 
running average for the data of a, b, and c. The arrows in Fig. [4] 
indicate the direction as the paths evolve with time. To follow 
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Table 3. Proper motions and apparent speeds for individual components (for further explanation see section 3.4). 



V 


Comp. 




/J ap p (coresep.) 




/3 app (rect. coord.) 


[GHz] 




[mas/year] 


M 


[mas/year] 


M 


5.0 


a 


-0.009±0.042 


-0.19+0.84 


0.117+0.058 


2.19+1.09 




b 


0.034±0.048 


0.68+0.95 


0.110+0.055 


2.06+1.03 




c 


0.045±0.017 


0.89+0.34 


0.121+0.060 


2.27+1.13 




d 


0.018±0.085 


0.36+1.68 


0.424+0.212 


7.96+3.98 


8.4 


a 


-0.025±0.018 


-0.49+0.37 


0.169+0.084 


3.17+1.58 




b 


-0.042±0.030 


-0.84+0.59 


0.205+0.102 


3.85+1.92 




c 


-0.014±0.046 


-0.29+0.91 


0.164+0.082 


3.08+1.54 




d 


0.016±0.031 


0.32+0.61 


0.276+0.138 


5.18+2.59 


15.3 


a 


-0.004±0.006 


-0.09+0.13 


0.24+0.12 


4.5+2.3 




b 


-0.011±0.009 


-0.22+0.19 


0.32+0.16 


6.0+3.0 




c 


-0.037±0.013 


-0.74+0.27 


0.54+0.27 


10.0+5.0 




d 


-0.029±0.037 


-0.58+0.73 


1.60+0.80 


29.8+15.0 


22.2 


a 


0.001 ±0.002 


0.02+0.05 


0.043+0.021 


0.80+0.40 




b 


-0.005±0.006 


-0.10+0.12 


0.088+0.044 


1.65+0.82 




c 


0.006±0.010 


0.14+0.20 


0.269+0.134 


0.50+0.25 




e 


0.019±0.006 


0.38+0.13 


0.005+0.002 


0.09+0.04 


43.2 


a 


0.015+0.012 


0.31+0.24 


0.031+0.015 


0.58+0.29 




b 


-0.011+0.014 


-0.23+0.27 


0.044+0.022 


0.82+0.41 




c 


-0.008+0.027 


-0.17+0.54 


0.179+0.089 


3.36+1.68 




d 


-0.016+0.021 


-0.32+0.43 


0.235+0.117 


4.41+2.20 




e 


-0.009+0.042 


-0.18+0.84 


0.775+0.387 


14.55+7.27 



the position angle evolution with time completely, more densely 
sampled data in time are required. 




I i i i i i i i I 

-0.6 -0.4 -0.2 

x 

Fig. 4. Motion of the inner four components at 15 GHz in rectan- 
gular coordinates (5-point running average for a, 4-point running 
average for b and c, and 3 -point running average for d). The ar- 
rows indicate the direction of the motion. 



3.4. Apparent velocities 

In the subsections above we showed that a "stationary" scenario 
- components showing no long-term outward motion - describes 
our results best. We calculate the apparent speeds based on the 
observed "stationarity" and list them in Table [3] The apparent 
speeds have been calculated for the inner 2.5 mas at 5.0 GHz, 
for the inner 3.5 mas at 8.4 GHz, for the inner 2.5 mas at 15.3 
GHz, for the inner 2 mas at 22.2 GHz, and for the inner 1.3 mas 
at 43.3 GHz. 

Some of the model components reveal apparent inward motion, 



which is not uncommon in AGN. We discuss this in the 
Discussion section. However, none of these values is of 3<x 
significance. For the apparent velocities derived for the motion 
based on the core separation evolution with time only, we 
thus exclusively find subluminal values. We list these values 
in Table [3] on the left side. The values have been calculated 
depending on the assumed angle to the line of sight of 4.9° or 
^0°, respectively. Angles to the line of sight are taken from 
Bach et al. (2005). 

The fast apparent velocities found by e.g., Bach et al. (see 
Table |4] taken from Bach et al. 2005) are not confirmed by the 
apparent speeds we derive for the "stationary" components. 
The MOJAVE 2cm survey webpage 

(http://www.physics.purdue.edu/astro/MOJAVE/sourcepages/) 
quotes a number of 1.199+ 0.116 mas/y for the motion in S5 
0716+714 (Lister et al., in prep.). This value is significantly 
higher than the values we find based on the core separation/time 
relation for individual jet components. 

In columns (5) and (6) of Table [3] we calculate lower limits 
for the proper motion and apparent speed based on the motion 
perpendicular to the radial motion - taking the motion with 
regard to the position angle into account. We estimated the 
distance traveled by each component in rectangular coordinates. 



3.5. Frequency dependence 

After having identified all components independently according 
to core separation, position angle, flux-density, and major axis 
- as described in previous sections - we then traced the indi- 
vidual components across the frequencies. We did not a priori 
make any assumptions concerning frequency dependent proper- 
ties. We averaged the component positions (weighted by their 
flux-density) and found that they can easily be grouped into the 
four components shown in Fig.|5ja). We there show the average 
core separation of each component at each frequency. The fre- 
quency identification is hampered by the different resolution at 
different frequencies. Thus, at 22 GHz, components close to the 
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Fig. 5. Identification of the components across the frequencies, 
(a) shows the average core separation of each component at each 
frequency, (b) shows the logarithm of the flux-density as func- 
tion of frequency, (c) shows the position of the components ac- 
cording to the frequency at which they were observed in rectan- 
gular coordinates along the jet. 



tially the right side of the jet and towards higher frequencies the 
left side of the jet becomes visible. 



4. Correlation: total flux-density 
evolution/position-angle evolution 

In Fig. 0a) we show the flux-density light-curves at 4.8, 8.0, 
and 14.5 GHz as obtained within the Michigan monitoring pro- 
gramme. Superimposed we show part of the epochs of the VLBI 
observations. This figure illustrates that morphological informa- 
tion is available to investigate a possible correlation between the 
longterm flux-density evolution and the kinematics. Although 
much more flux-density compared to VLBI information is avail- 
able, the VLBI information covers the general trend well. In par- 
ticular, maxima as well as minima in the flux-density light-curve 
can be accompanied by interferometric measurements. 
In Fig. |6jb) we show a collection of several data sets (all at 15 
GHz): the total flux density at 14.5 GHz (black dots, data), the 
flux density of the core (dot-dashed), the flux density of compo- 
nent a (closest component to the core, dashed), and the position 
angle for component a (dotted) as function of time. We find that 
most of the total flux density is contained in the core compo- 
nent - which is usual for this type of object. We also find evi- 
dence for a correlation between the total flux-density evolution 
with time and the position-angle evolution of component a with 
time. Thus, we find a correlation between flux-density variabil- 
ity and kinematic properties for the innermost jet component. 
This means as well that the "core" is exhibiting behavior similar 
to a jet-component. This might suggest that the longterm flux- 
density evolution in S5 0716+714 may have a non-negligible 
geometric contribution. The position angle information for a is 
not sufficient to unambiguously prove a correlation, but the data 
consistently reproduce the same trend as seen in the total flux- 
density light-curve. We plan to probe and investigate further this 
particular relation in more detail with more data and on shorter 
timescales. 



5. Discussion 

5.1. "Fast" versus "Slow" scenario 

In this paper we show that motion in the pc-scale jet of 0716+714 
can be interpreted within a new kinematic scenario of non-radial 
motion. Components do not show the generally assumed fast ap- 
parent outward motion. Instead, stationary components with re- 
gard to the core separation are observed. In addition, significant 
motion perpendicular to the jet is observed. In the following, we 
discuss our results and compare them to those from other studies 
of AGN pc-scale jet kinematics. 



core can be seen, which are not visible at lower frequencies. We 
want to mention that - for the afore mentioned reason - at 15 GHz 
we grouped the components a and b into one component and at 
22 GHz the components c and d into one component and e and 
xi and X2 and X3 into one component. We find that we can easily 
identify four components (labeled I-IV) between 5 and 15 GHz. 
In Fig.|3b) we show the logarithm of the flux-density as a func- 
tion of frequency for these 4 components. Fig. [5jc) shows the 
position of the components according to the frequency at which 
they were observed in rectangular coordinates along the jet. We 
can trace components at different frequencies along the jet and 
find that at lower frequencies we see consistently and preferren- 



5.2. Apparent stationarity and non-radial motion of jet 
components 

Our observations reveal a new picture of jet component motion 
in the BL Lac Object 0716+714: We do not find any longterm 
outward motion at any of the investigated frequencies. Instead, 
model components remain at similar core separations in the in- 
vestigated time span. This is in contrast to previous studies of 
jet component motion in AGN, where components tend to show 
predominantly outward motion (e.g., Britzen et al. 2007, 2008). 
In contrast to previous studies — that relied mainly on investi- 
gating the core separation evolution — we study the position- 
angle evolution with time as well. Although we do not find any 
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evidence for long term outward motion, we instead find that all 
modelled jet components move with regard to the position an- 
gle. A similar motion phenomenon has already been observed 
and analyzed in the jet of 1803+784 (Kudryavtseva et al. 2007, 
Britzen et al. 2009). 

In quite a number of AGN, jet components follow either the 
same path or a small number of different paths, which have been 
explained within precession scenarios (e.g., PKS 0528+134: 
Britzen et al. 1999). While components in general tend to fol- 
low the same jet ridge line in a number of AGN investigated in 
detail, the jet ridge line itself in 0716+714 shows strong evolu- 
tion with time. 

This can be explained as a result of precession of the base 
of the jet. This precession can be explained and described by 
a number of different models. Assuming that the core of the 
AGN hosts a single object (i.e. SMBH), precession can oc- 
cur as a result of the Bardeen-Petterson effect (Caproni et al. 
2006), magnetic torques (Lai 2003) or magnetohydrodynamic 
instabilities, prominently Kelvin-Helmholtz instabilities (e.g., 
Camenzind & Krockenberger 1992, Hardee & Norman 1988, 
Birkinshaw 1991, Zhao et al. 1992, Hardee et al. 1994, Hardee 
et al. 1997, Meier & Nakamura 2006, Perucho et al. 2006). 
If we assume that the AGN core hosts a binary system (presum- 
ably a binary black hole, BBH) then precession can occur as a 
result of the orbital motion of the binary system (e.g., Lobanov 
& Roland 2005, Roos et al. 1993, Kaastra & Roos 1992), tidal 
forces, gravitational torques (e.g., Katz 1997, Romero et al. 
2000), as well as orbital motion of the system around the galac- 
tic gravity center (Roland et al. 2008). 

Finally, it is possible that the bending of the jet ridge line is 
explained by external causes. Interaction with the surrounding 
intergalactic medium can cause a bending of the jet, it is, how- 
ever, highly unlikely that this bending will show any periodic 
or quasi-periodic behaviour. Furthermore, bending of a jet ridge 
line can occur by gravitational interaction of the host galaxy 
with a nearby galaxy (Blandford & Icke 1978, Jagers & de Grijp 
1985, Lupton & Gott 1982). Again, each and all the above phe- 
nomena, of external origin, can only partially but not fully ex- 
plain the bending of the jet ridge line. For a more detailed dis- 
cussion of the above models please see Britzen et al. (2009) and 
references therein. 

We are currently working on a model to explain the quasi- 
stationarity within a model of non-ballistic motion (Gong 2008). 



5.3. TeV blazars and subluminal motion 

In quite a large fraction of the TeV blazars, subluminal motion 
has been detected (e.g., Piner et al. 2008). The apparent pattern 
speeds measured in the TeV blazars are considerably slower 
than those measured in sources selected for their compact radio 
emission (e.g., CJF: Britzen et al. 2008), or for their GeV 
gamma-ray emission (Jorstad et al. 2001). S5 0716+714 has 
most recently been detected in the TeV-regime as well. The 
subluminal motion we find in this source thus fits nicely with 
the tendency for significantly slower apparent speeds in this 
class of objects. Piner et al. (2008) discuss the possible origin 
of these slow apparent motions and their connection to the 
TeV production mechanism. Since the rapid radio variability 
detected in S5 0716+714 requires quite high Doppler factors 
(e.g., Wagner & Witzel 1995 and references therein), it is likely 
that the slow apparent speeds we see result from pattern speeds 
and do not reflect the speeds of the underlying bulk motion. 



5.4. Correlation between kinematics & flux-density 
variability 

We find some evidence for a possible correlation between po- 
sition angle changes of the innermost component and the long- 
term flux-density variability. This favors a geometric contribu- 
tion to the origin of the longterm flux-density evolution. Further 
support for the geometric contribution has been discussed based 
on monitoring of the flux-density variability and studies of the 
light-curves: 

The flux-density variability in the optical and in the radio band 
shows different timescales (Raiteri et al. 2003). Raiteri et al. find 
that the variation of the amplitude of the flux-density decreases 
with increasing wavelength and that the outbursts in the optical 
and radio regime are caused by the same mechanism. The flux- 
density in the radio regime seems to be the sum of two different 
spectral distributions. Raiteri et al. claim that the oscillations in 
both wavelength regimes have either an energetic or a geometric 
origin. Nesci et al. (2005) find in an optical light curve some ev- 
idence for long term variability on which short term variations 
are superimposed. Indications for a precessing jet have also been 
found (Nesci et al.). 

The tentative correlation that we find between the position an- 
gle changes of the innermost component and the long-term flux- 
density variability leads us to believe that a geometrical expla- 
nation for the long-term variability of S5 0716+714 is plausi- 
ble. Camenzind & Krochenberger (1992), first introduced the 
lighthouse effect, according to which the periodiocity (or quasi- 
periodicity) observed in the lightcurves of blazars can be at- 
tributed to the fast rotation of non-axisymmetric plasma bub- 
bles near the central black hole. During the propagation of the 
bubble downstream in the jet, it experiences beaming effects, as 
the viewing angle changes, thus producing the outbursts in the 
light curves. The possible quasi-periodicity found from the light 
curves is characteristic of the rotation period of the plasma near 
the black hole, although it is significantly shorter due to rela- 
tivistic effects. The above scenario requires a well collimated jet 
with a relatively small opening angle. In any different case the 
flux-density dies out quickly. 

Schramm et al. (1993), find that the lighthouse effect is possi- 
bly responsible for the optical variability for the blazar 3C 345, 
although they consider that extrinsic effects e.g. microlensing) 
could also contribute to the long-term variability observed. For 
OJ 287, Katz (1997) argues that a precessing disk consequently 
results in beaming effects for the jet and variability in the flux- 
density. For the same source, Villata et al. (1998), proposed a 
beaming model for which both black holes of a BBH system 
produce a jet. A combination of bending and beaming of the 
jets can produce the double-peaked lightcurve observed for that 
source, thus attributing the variability only to geometrical ef- 
fects. Qian et al. (2007) make use of a similar model to explain 
the light curve of quazar 3C 345.3. Rieger (2004) also attributes 
variability of blazars to differential Doppler boosting effects, as 
the plasma follows helical trajectories, induced by different pos- 
sible mechanisms (as orbital motion of a BBH system, intrinsic 
jet rotation and jet precession). A geometrical origin (relativistic 
beaming) for the variability of blazar AO 0235+16 is also con- 
sidered by Raiteri et al. (2001). 

Non geometrical phenomena can, however, also explain the 
long-term variability seen in the light curves of 0716+714. Most 
prominently, two different models have been proposed for OJ 
287 which shows periodic variability in its light curves (in both 
the radio and the visual regime). Sillanpaa. et al. (1988) orig- 
inally argued that tidal forces from the secondary black hole 
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Fig. 6. (a) shows the flux-density light curves obtained by the Michigan monitoring programme at three different frequencies. 
Superimposed are part of the epochs of the available VLBI epochs investigated in this paper. In (b) a compilation of several datasets 
is displayed: the total flux-density at 14.5 GHz (black dots with error bars), the flux density of the core (dot-dashed), the flux 
density of component a (dashed) as function of time and the position angle of component a (dotted). Most of the total flux-density 
is contained in the core component, as expected. 



around the pericenter of its orbit will lead to an increased accre- 
tion rate for the primary black hole and thus result in an outburst. 
Lehto & Valtonen (1996) proposed that the secondary black hole 
plunges into and through the accretion disk of the primary before 
and after the pericenter, resulting in the double outburst observed 
in the light curves (see also Valtonen et al. 2006). This proce- 
dure however, as argued by Valtaoja et al. (2000), is thermal and 
should not produce an outburst in the radio regime; thus it seems 
unlikely for the case of 0716+714. In summary, we find that a 
geometrical explanation is a better match to the observational 
properties of 0716+714. 



6. Conclusions 

In the previous sections we presented our results on the kine- 
matics of jet component motion in S5 0716+714. We find a 
different scenario than previously described in the literature for 
this source. In particular, we find no long term outward motion 
for this source. We present a new model component motion 
scenario for S5 0716+714 with the following implications: 

• The jet ridge line changes significantly from epoch to epoch. 
Components do not follow a well-defined path. 

• In this scenario the jet components show no 3cr evidence for 
long term outward motion. Instead, they remain at similar core 
separations and seem to change significantly in position angle 
within the computed uncertainties. We calculate exclusively 
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apparent subluminal speeds based on the observed core separa- 
tion/time relation. 

• S5 0716+714 has been detected as a TeV blazar and shows 
- as most of the other TeV blazars - subluminal motion. The 
observed apparent motions most likely result from pattern 
speeds. 

• Taking the "real" motion in rectangular coordinates into 
account, we find significantly higher values for the apparent 
speeds. These values can only be limits. 

• We find a possible correlation between the total flux density 
and the position angle of the component closest to the core 
(a). This relation might suggest a non-negligible geometric 
contribution to the origin of the longterm flux-density variability 
in S5 0716+714. 

The usually observed outward motion of jet components in 
AGN is — according to our analysis — not observable in 
the BL Lac object S5 0716+714. Instead, we find kinematic 
properties, e.g. non-radial motion, evolving jet ridge-line, etc. 
which have not been reported for this source so far and have 
only been observed in a couple of other AGN so far (e.g. Britzen 
et al. 2009). More observations and modelling are required to 
understand the physical origin of these phenomena. 

Acknowledgements. We highly appreciate many comments and suggestions 
by T.P. Krichbaum on how to check the significance of the results presented 
here. I. Agudo acknowledges financial support from the EU Commission 
for Science and Research under contract HPRN-CT-2002-00321 (ENIGMA 
Network). M. Karouzos was supported for this research through a stipend from 
the International Max Planck Research School (IMPRS) for Radio and Infrared 
Astronomy. We wish to acknowledge the 2cm Survey/MOJAVE team. This re- 
search has made use of data from the University of Michigan Radio Astronomy 
Observatory which has been supported by the University of Michigan and the 
National Science Foundation. This work made use of the VLBA, which is a fa- 
cility of the National Science Foundation, operated under cooperative agreement 
by Associated Universities, Inc. We used the European VLBI Network, which is 
a joint facility of European, Chinese, South African and other radio astronomy 
institutes funded by their national research councils. This work is also based on 
observations with the 100m radio telescope of the MPIfR (Max-Planck-Institut 
fur Radioastronomie) at Effelsberg. 



References 

Agudo, I, Bach, U., Krichbaum, T. P., Marscher, A. P., et al. 2007, A&A, 476, 
17 

Agudo, I., Krichbaum, T.P, Ungerechts, H., Kraus, A., et al. 2006, A&A, 456, 
117 

Anderhub, H., Antonelli, L.A., Antoranz, P., Backes, M., et al. 2009, ApJ, 704L, 
129 A 

Antonucci, R.R.J., Hickson, P., Olszewski, E.W., Miller, J.S. 1986, AJ, 92, 1 
Bach, U., Krichbaum, T.P, Ros, E., et al. 2005, A&A, 433, 815 
Birkinshaw, M. 1991, MNRAS, 252, 505 
Blandford, R.D. & Icke, V. 1978, MNRAS, 185, 527 
Britzen, S., Witzel, A., Krichbaum, T.P, et al. 2005, MNRAS, 362, 966 
Britzen, S., Vermeulen, R.C., Taylor, G.B., Campbell, R.M., et al. 2007, A&A, 
472, 763 

Britzen, S., Vermeulen, R.C., Campbell, R.M., et al. 2008, A&A, 484, 119 

Britzen, S., Kudryavtseva, N, Witzel, A., et al. 2009, A&A, accepted 

Britzen, S., Meyer, V., Witzel, A., et al., Proceedings of the 8th European VLBI 

Network Symposium. September 26-29, 2006,Torun, Poland 
Camenzind, M. & Krockenberger, M. 1992, A&A, 255, 59 
Caproni, A., Abraham, Z., Mosquera Cuesta, H.J. 2006, ApJ, 638, 120 
Eckart, A., Witzel, A., Biermann, P., et al. 1987, A&AS, 67, 121 
Fey, A.L., & Chariot, P. 2000, ApJS, 128, 17 

Gabuzda, D.C., Kovalev, Y.Y., Krichbaum, T.P, et al. 1998, A&A, 333, 445 

Gong, B.P 2008, MNRAS, 389, 315 

Hardee, P.E., Norman, M.L. 1988, ApJ, 334, 70 

Hardee, P.E., Cooper, M.A., Clarke, D.A. 1994, ApJ, 424, 126 

Hardee, P.E., Clarke, D.A., Rosen, A. 1997, ApJ, 485, 533 

Jagers, W.J. & de Grijp, M.H.K. 1985, A&A, 143, 176 

Jorstad, S.G., Marscher, A.P., Mattox, J.R., et al. 2001, ApJS, 134, 181 

Kaastra, J.S. & Roos, N. 1992, A&A, 254, 96 



Katz, J.I. 1997, ApJ, 478, 527 

Kellermann, K.I., Vermeulen, R.C., Zensus, J.A., & Cohen, M.H. 1998, AJ, 1 15, 
1295 

Kellermann, K.I., Lister, M.L., Homan, D.C., et al. 2004, ApJ, 609, 539 
Kudryavtseva, N, et al., Proceedings of the 8th European VLBI Network 

Symposium. September 26-29, 2006, Torun, Poland 
Lai D. 2003, ApJ, 591, 119 
Lehto, H.J., Valtonen, M.J. 1996, ApJ, 460, 207 
Lobanov, A.P & Roland, J. 2005, A&A, 431, 831 
Lister, M., & Homan, D. 2005, AJ, 130, 1418 

Lister, M.L., Cohen, M.H., Homan, D.C., Kadler, M. et al. 2009, AJ, accepted, 
arXiv0909.5100 

Lovell, A.J. 2000, Astrophysical Phenomena Revealed by Space VLBI, 
Proceedings of the VSOP Symposium, H. Hirabayashi, PG. Edwards, and 
D.W. Murphy (eds.), Institute of Space and Astronautical Science, 301 
Lupton, R.H., Gott, J.R., III 1982, ApJ, 255, 408 
Meier, D.L., Nakamura, M. 2006, ASP Conference Series, Vol. 350 
Meyer, V.A. 2007, Diploma Thesis, University of Cologne 
Nesci, R., Massaro, E., Rossi, C, Sclavi, S., et al. 2005, AJ, 130, 1466 
Nilsson, K., Pursimo, T, Sillanpaa, A., Takalo, L.O., Lindfors, E. 2008, A&A, 
487, 29 

Perez-Torres, M.„ Marcaide, J., Guirado, J., et al. 2004, A&A, 428, 847 
Perucho, M., Lobanov, A.P, Marti, J.M., Hardee, PE. 2006, A&A, 456, 493 
Piner, G.B., Pant, N, Edwards, P.G. 2008, ApJ, 678, 64 
Quirrenbach, A., Witzel, A., Wagner, S., et al., 1991 ApJ, 372, 71 
Raiteri, CM., Villata, M., Aller, H.D., Aller, M.F., et al. 2001, A&A, 377, 396 
Raiteri, CM., Villata, M., Tosti, G, et al., 2003 A&A, 402, 151 
Rastorgueva, E.A., Wiik, K., Savolainen, T., Takalo, L.O., et al. 2009, A&A, 
494, 5 

Rieger, FM. 2004, ApJ, 615, L5 

Roland, J., Britzen, S., Kudryavtseva, N.A., Witzel, A., Karouzos, M. 2008, 
A&A, 483, 125 

Romero, C.E., Chajet, L., Abraham, Z., Fan, J.H. 2000, A&A,360, 57 
Roos, N., Kaastra, J.S., Hummel, C.A. 1993, ApJ, 409, 130 
Ros, E., Marcaide, J.M., Guirado, J.C, & Perez-Torres, M.A., 2001 A&A, 376, 
1090 

Schalinski, C.J., Witzel, A., Krichbaum, TP, et al. 1992, in proceedings of 
Variability of Blazars, E. Valtaoja & M. Valtonen (eds.) (Cambridge CUP), 
225 

Schramm, K.-J., Borgeest, U., Camenzind, M., Wagner, S.J., et al. 1993, A&A, 
278, 391 

Sillanpaa, A., Haarala, S., Valtonen, M.J., Sundelius, B., et al. 1988, ApJ, 325, 
628 

Teshima, M., The MAGIC Collaboration 2008, ATel 1500, IT 

Valtaoja, E., Tersranta, H, Tornikoski, M., Sillanpaa, A. 2000, ApJ, 531, 744 

Valtonen, M.J., Lehto, H.J., Sillanpaa, A., Nilsson, K, et al. 2006, ApJ, 646, 36 

Villata, M., Raiteri, CM., Sillanpaa, A., Takalo, L.O. 1998, MNRAS, 293, 13 

Wagner, S.J. & Witzel, A. 1995, ARA&A, 33, 163 

Witzel, A., Schalinski, C.J., Johnston, K.J., et al. 1988, A&A, 206, 245 

Zhao, J.-H, Burns, J.O.; Hardee, P.E., Norman, M.L. 1992, ApJ, 387, 69 

Zensus, J.A., Ros, E., Kellermann, K.I., et al. 2002, AJ, 124, 662 



S. Britzen: Non-radial motion in the TeV blazar S5 0716+714 



13 




Fig. 7. 5 GHz VLBI images of S5 0716+714. Model-fits convolved with the interfero metric beam and with the residual map added 
to them are shown with the Gaussian model components superimposed. The following figures show similar images obtained at 
different frequencies. 
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Fig. 8. Images of S5 0716+714 at 8 GHz. 
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Fig. 9. First part of images of S5 0716+714 at 15 GHz. 
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Fig. 10. Second part of images of S5 0716+714 at 15 GHz. 
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Fig. 11. First part of images of S5 0716+714 at 22 GHz. 




Fig. 12. Second part of images of S5 0716+714 at 22 GHz. 



S. Britzen: Non-radial motion in the TeV blazar S5 0716+714 



19 




Fig. 13. Images of S5 0716+714 at 43 GHz. 



20 S. Britzen: Non-radial motion in the TeV blazar S5 0716+714 

Table 5. The parameters of the model-fitting. 



v Epoch 


S 


r 


e 


M.A. Comp. 


[GHz] 


[mJy] 


[mas] 


[cleg] 


[mas] 



5.0 1992. 
1992. 
1992. 
1992. 
1992. 



73 
73 
73 
73 
73 



600+30 
35+2 
12+2 
10+2 
5+1 





0.69+0.03 
1.98+0.10 
4.25+0.21 
6.52±0.33 



0. 

24.1 + 1.20. 
7.1 ±0.40. 
18.0+1.00. 
12.0+0.6 0. 



14+0.01 
01±0.01 
01±0.01 
53+0.03 
51+0.03 



Core 
a 
c 
e 
f 



1994.70 78+4 0.01±0.01 Core 

1994.70 231 + 12 0.32±0.02 1.8+0.1 0.22±0.01 a 
1994.70 14+2 1.87+0.10 15.5+0.1 0.77±0.04 c 
1994.70 16+3 3.32±0.20 12.9+0.6 0.89±0.04 d 
1996.63 182+10 0.03+0.01 Core 

1996.63 22+1 1.15+0.06 6.5±0.3 0.01+0.01 b 
1996.63 5±1 3.75+0.2020.2+1.00.25±0.01 d 

1996.63 7.0+0.4 5.66±0.30 23.0+1.2 1.78±0.09 f 
1999.89 552+30 0.15+0.01 Core 

1999.89 40+4 1.12+0.06 8.6+0.40.01+0.01 b 
1999.89 20+4 2.30±0. 11 12.3+0.6 0.62±0.03 c 
1999.89 12+2 3.98±0.20 12.8+0.6 1.70±0.09 d 

Core 
a 
b 
c 
d 
e 
f 



2000. 
2000. 
2000. 
2000. 
2000. 
2000. 
2000. 



82 
82 
82 
82 
82 
82 
82 



446+22 
48+2 
25+2 
12+2 
7.0+0.4 
4.0+0.2 
3±1 





0.51 ±0.03 
1.36+0.07 
2.23+0.11 
3.25+0.16 
4.32+0.22 
6.18+0.31 



0. 

26.9+1.3 0. 
13.8+0.7 0. 
11.7+0.60. 
11.4+0.60. 
11.6+0.60. 
16.9+0.8 2. 



07+0.01 
01±0.01 
29±0.01 
52+0.03 
85+0.04 
97+0.05 
17+0.11 



8.4 1994.21 261 + 14 0.14+0.01 Core 

1994.21 26+13 0.79+0.0413.8+0.7 0.45+0.02 a 
1994.21 10+15 1.82±0.09 16.7+0.8 0.75+0.04 b 
1994.21 12+24 3. 31+0.17 12.8+0.6 1.15±0.06 d 
1995.65 238+12 0.07+0.01 Core 

1995.65 47.00+2.35 0.370±0.02 9.9+0.5 0.23+0.01 a 
1995.65 7.0+0.7 1.37+0.0715.0+0.8 0.49+0.02 b 
1995.65 3.0+0.6 2.63±0.13 19.4+1.00.36+0.02 c 
1995.65 4.00±0.08 3. 86±0.19 15.5+0.8 1.05±0.05 d 
1997.03 191 + 10 0.05+0.01 Core 

1997.0321.00±0.11 0.52±0.03 22.5+1. 1 0.01±0.01 a 
1997.93 372+19 0.02+0.01 Core 

1997.93 31.0+1.6 0.54+0.03 10.7±0.5 0.36+0.02 a 
1997.93 15.0+1.5 1.60+0.08 8.5+0.4 0.34+0.02 b 
1997.93 7.0+1.4 3.24±0.16 15.8+0.8 0.86±0.04 d 

Core 
a 
b 
c 
d 



1999. 
1999 
1999 
1999 
1999 



41728.0+36.4 
41233.0+11.7 
41 32+16 
41 11.00+0.22 
41 7.0+1.4 





0.24+0.01 
0.99+0.05 
2.05+0.10 
3.61+0.18 



0. 

13.2+0.7 0. 
13.6+0.7 0. 
9.0+0.5 0. 
12.6+0.6 0. 



01±0.01 
12+0.01 
30+0.02 
45+0.02 
90±0.05 



2002.48 


663+33 


0.01+0.01 


Core 


2002.48 


152+1 


0.24+0.01 27.5+0.20.15+0.01 


a 


2002.48 


17.0+0.5 


1.15+0.01 18.2+0.10.46+0.02 


b 


2002.48 


7.0+0.7 


2.20+0.01 19.7+0.20.52+0.03 


c 


2002.48 


10+2 


3.69±0.18 13.6+0.1 1.08±0.05 


d 



2003. 
2003 
2003 
2003 
2003 



20 
20 
20 
20 
20 



1106±60 
31+2 
23+2 
4.0+0.4 
4.0+0.8 





0.58+0.03 
1.55+0.08 
2.67+0.13 
3.58+0.18 



0. 

32.1 + 1.60. 
19.7+1.00. 
18.2+0.9 0. 
16.5+0.8 0. 



12±0.01 
01±0.01 
76+0.04 
01±0.01 
70±0.04 



Core 
a 
b 
c 
d 
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v Epoch 


c 


r 9 M.A. Comp. 


[GHz] 


[mjy] 


[mas] [deg] [mas] 




2003 20 


4.0+0.5 


A ^9-lO 09 99 O-i-O Q O /19-lO 09 
4.JZ+U.UZ ZZ. \}±\}.y U.4Z+U.UZ 


e 


2003.49 


1518+80 


O O O 09_i_0 01 

U UU.Uz+U.Ul 


core 


2003.49 


96+14 


o ^7x0 O/l /I _i_o in ni _i_r\ m 
U.J /+U.U4 JJ.4+U.J U.U1+U.U1 


a 


2003.49 


11 + 1 


1 1 ^_i_0 f\< 90 Q_i_0 /I O 01 _i_0 01 

1 . 1 j+U.Uj zU.y+U.4 U.U1 +U.U1 


u 
u 


2003.49 


16+1 


9 30+0 19 91 4+0 4 1 90+0 06 




2003.49 


6±2 


t-. jjiU.UD Ld.L± 1 .Z Z.UO±U. 1U 


e 


2003.49 


3.0+0.5 




X 


15 3 1994 67 


321+20 


0.05+0.01 


Core 


1994.67 


71+4 


0.17±0.01 4.6+0.20.01+0.01 


a 


1994.67 


20+1 


0.80±0.04 34.0+ 1 .7 0.84+0.04 


b 


1994.67 


9.0+0.9 


1.57±0.08 11.9+0.6 0.58±0.03 


c 


1996.53 


208+10 


0.05+0.01 


Core 


1996.53 


26.0+1.3 


0.23+0.01 17.6+0.9 0.01+0.01 


a 


1996.53 


19+10.67±0.034 6.2±0.3 0.28±0.01 


b 


1996.53 


j.U+U.j 


1 ztQj-0 OR 1 9 £-i-0 £ O 1 1 4-0 09 


c 


1996.53 


i o+o 9 


9 41 -i-O 19 18 £-1-0 Q O 48-1-0 09 
Z.+ 1+U.1Z lo.O±U.7 U.^ru+U.UZ 


d 


1996.53 


2.0+0.4 


4 1 5+0 21 16 3+080 67+0 04 


X 


1996.82 




01 -i-O 01 


Core 


1996.82 


1 ^+1 
1 J±l 


^8+0 09 1 7 ^+0 Q 01 -i-O 01 


a 


1996.82 


8.0+0.8 


1 48+0 01 1 1 0+0 40 53+0 03 


c 


1999.55 


i o^+^9 


07-1-0 01 


Core 


1999.55 


1 lA-t-Q 
1 / ++y 


O 1 8-1-0 01 97 9-1-1 AO 01 4-0 01 
U.loxU.Ul Z/.Z+1.4 U.UlxU.Ul 


a 


1999.55 


99-1-1 
ZZ± 1 


76-1-0 04 1 ^ 6-1-0 7 OQ-i-0 01 


b 


1999.55 


1 6 0+ 1 6 


1 SS+O 08 1 1 S+O 6 ^6+0 09 


c 


1999.55 


5+1 


3 24+0 16 17 8+090 36+0 02 


e 


2001.17 


^4s+^o 

J+J± JU 


OS+O 01 


Core 


2001.17 


79 j-99 
/ Z+ZZ 


O zt9j-0 09 91 ^-i-1 OO 1 6-1-0 01 
U.^fZ+U.UZ Zl.J+l.UU. 10+U.vl 


a 


2001.17 


9 0-1-0 ^ 


8Q-I-0 09 1 6 8-1-0 ^0 01 -i-O 01 


b 


2001.17 


19+^ 


1 84-1-0 01 1 7 O-i-O 10 61 -i-O 0^ 
l.OTxW.Ul 1 / .UxU.l U.UlxU.Uj 


c 


2001.17 


^ 0-1-1 ^ 
J.U+ 1 .J 


^ 1 O-i-O 01 1 ^ ^4-0 7 6Q4-O O'* 


e 


2001.17 


9 O-i-O ^ 


S 90-1-0 09 1 7 4-1-0 Q 01 -i-O 01 
J.Zuxu.UZ 1 / .+±U .y W .W 1 ±W.W1 


X 


2003 65 


99QQ+1 i o 


09+0 01 

\J \J \J .\J ZjZL.\J .\J L 


Core 


2003 65 


254+13 


O 16+0 0] ,6 0+1 8 10+0 01 

\J a J. \J -!—\J .Ul *J\J .\J _!_ -L a U \J . ±\J -!—\J .\J -I 


a 


2003.65 


15.0+0.8 


0.70+0.04 32.7+1 .6 0.33+0.02 


b 


2003.65 


10+2 


2.20+0.11 18.1+0.91.03+0.05 


d 


2004.44 


1403±70 


00.01±0.01 


Core 


2004.44 


419±20 


0.17±0.01 29.2+1.5 0.07+0.01 


a 


2004.44 


84.0+4.2 


0.56+0.03 31. 1 + 1.60.18+0.01 


b 


2004.44 


25+3 


1.03±0.05 24.2+1.2 0.25±0.01 


c 


2004 60 

_ w v./ r . v./ v./ 


1336±70 


00.01±0.01 


Core 


2004 60 

z— v./ v./ r 


200+10 


0.26+0.01 25.5+1.3 0.01±0.01 


a 


2004.60 


13+7 


1.30±0.07 29.9+1.5 0.31 ±0.02 


c 


2004.79 


1841 ±90 


0.07±0.04 


Core 


2004.79 


15.0+0.8 


0.33+0.02 20.8+1.00.01+0.01 


a 


2004.79 


100+5 


0.69±0.04 21.4+1 .0 0.44+0.02 


b 


2004.79 


8.0+0.8 


1.46+0.07 17.8+0.9 0.31+0.02 


c 


2004.79 


3.0+0.5 


2.02±0.1021.7+1.1 0.81±0.04 


d 


2004.97 


1018+43 


0.01±0.02 


Core 


2004.97 


288+40 


0.16+0.01 24.6+0.20.07+0.01 


a 


2004.97 


51.0+0.5 


0.80±0.04 18.6+1.00.32+0.02 


b 


2004.97 


17+3 


1.31±0.0421.9+0.3 0.09+0.01 


c 


2004.97 


2.0+0.5 


1.94+0.03 22.8+0.40.01+0.01 


d 


2004.97 


3.0+0.2 


3.28+0.08 33.4+1.7 0.01+0.01 


e 


2004.97 


2+1 


4.16±0.02 22.8+0.1 0.01±0.01 


X 


2005.01 1406.0+0.5 


0.05+0.01 


Core 
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v Epoch 


S 


r 6 M.A. Comp. 


[GHz] 


[mJy] 


[mas] [deg] [mas] 




2005.01 


91+4 


0.27+0.0420.4+0.1 0.01+0.01 


a 


2005.01 


16+8 


0.43+0.01 18.3+0.2 0.11+0.01 


b 


2005.01 46.00+0.05 


1.05+0.05 23.8+1.2 0.44+0.02 


c 


2005.01 


5.00±0.03 


1 , a r\n o h ■ a a a ai . a ai 

1.99±0.09 8.7+0.4 0.02+0.01 


A 

d 


2005.09 


319+16 


0.01+0.01 


Core 


2005.09 


84.0+4.2 


0.24+0.01 14.1+0.7 0.15+0.01 


a 


2005.09 


30.0+1.5 


0.86+0.0425.0+1.3 0.36+0.02 


b 


2005.09 


12.0+1.2 


1.46+0.07 19.2±1.0 0.10+0.01 


c 


2005.09 


12±2 


2.71+0.1416.9+0.8 1.65+0.08 


d 


2005.09 


2.0+0.6 


C /I O i A O OTiA/1 AA1 iAA1 

5.48+0.3 8.3+0.4 0.01+0.01 


X 


2005.70 


868+50 


0.03+0.01 


Core 


2005.70 


292+12 


0.12+0.0127.5+1.7 0.05+0.01 


a 


2005.70 


98.0+5.5 


0.57+0.03 18.9+0.8 0.22+0.01 


b 


2005.70 


9.0+22.5 


1 11 , ( \ rts- 1 A 1 i A A A 11 ■ A A1 

1.12+0.06 19.2+0.4 0.22+0.01 


c 


2005.70 


5+2 


2.00+0.29 19.0+0.9 0.68+0.03 


d 


2005.70 


5.0+1.5 


3.43+0.51 16.1+0.8 0.96+0.05 


e 


2005.72 


769+50 


0.01+0.01 


Core 


2005.72 


87.0+5.5 


0.57+0.03 20.5+0.8 0.21±0.01 


b 


2005.72 


14.0+2.5 


1.00+0.06 19.9±0.4 0.01+0.01 


c 


2005.72 


9+2 


2.58+0.2917.3+0.9 1.44+0.07 


A 

u 


2006.32498.0+24.9 


0.04+0.01 


core 


2006.32 


53.0±2.6 0.280+0.002 19.3+1.0 0.10+0.01 


a 


2006.32 


10.0±0.5 


0.71+0.0126.9+1.3 0.01+0.01 


b 


2006.32 


5.0+0.5 


1.26+0.0621.9+1.1 0.12+0.01 


c 


2006.32 


13.0+1.9 


1.85+0.0912.7+0.6 0.61+0.03 


A 

a 


2006.32 


2.0+0.4 


3.39±0.17 4.4+0.2 0.01+0.01 


e 


22.2 1992.85 


700+35 


0.13+0.01 


Core 


1992.85 


20+3 


1.62+0.08 7.2+0.4 0.01+0.01 


e 


1992.85 


10+2 


2.42+0.1217.9+0.9 0.01+0.01 


x2 


1993.71 


248.0±1.2 


0.01+0.01 


Core 


1993.71 


99+5 0.170+0.009 16.7+0.8 0.10+0.01 


a 


1993.71 


17.0+0.9 


0.57+0.03 15.3+0.8 0.13+0.01 


b 


1993.71 


16+16 


1.23+0.06 20.3+1.0 0.24+0.01 


d 


1994.21 


252.0±1.3 


0.01+0.01 


Core 


1994.21 


44+22 


0.18+0.01 22.1 + 1.1 0.01+0.01 


a 


1994.21 


13.0+0.7 


0.80+0.0418.2+0.9 0.37±0.02 


c 


1994.21 


6.0+0.9 


1.80+0.0915.2+0.8 0.47+0.02 


e 


1995.15 677.0+33.9 


0.05+0.01 


Core 


1995.15 


72+36 


0.17±0.01 4.1+0.2 0.18+0.01 


a 


1995.15 


4.00+0.04 


1.15+0.0613.6+0.7 0.01±0.01 


c 


1995.15 


5+1 


2.63+0.13 7.9+0.4 0.01+0.01 


x3 


1995.15 


8+16 


3.64±0.18 16.2+0.8 0.14+0.01 


y 


1995.47 


173+9 


0.04+0.01 


Core 


1995.47 


11.0+0.6 


0.54+0.03 17.2+0.9 0.01+0.01 


b 


1995.47 


7+14 


2.12+0.11 15.5+0.8 0.01+0.01 


xl 


1995.65 


174.0±8.7 


0.06+0.01 


Core 


1995.65 


123.0+6.2 


0.14+0.01 10.6+0.5 0.06+0.01 


a 


1995.65 


28+14 


0.36+0.0212.1+0.6 0.23+0.01 


b 


1995.65 


4.0+0.2 


0.76+0.04 10.2±0.5 0.01+0.01 


c 


1995.65 


5.0+0.8 


1.61±0.08 17.3+0.9 0.43+0.02 


e 


1996.34 


178+9 


0.01+0.01 


Core 


1996.34 


57.0+0.3 


0.18±0.01 10.5+0.5 0.01+0.01 


a 


1996.34 


24+12 


0.58+0.03 5.1+0.3 0.15+0.01 


b 


1996.34 


13.0+2.6 


1.99+0.1023.2+1.2 1.23 +0.06 


xl 


1996.90 


218+11 


0.01+0.01 


Core 
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v Epoch 


c 

o 


r 6 M.A. Comp. 


[GHz] 


[mJy] 


[mas] [deg] [mas] 




1996.90 


61.0+3.1 


0.14±0.01 21.5+1.1 0.13±0.01 


a 


1996.90 


7.00+0.35 


/ \ ill . r\ r\A lAO.rvCA r\C\ , A A1 

0.81±0.04 10.8±0.5 0.09±0.01 


c 


1996.90 


6.0+0.9 


1.57±0.08 10.9±0.5 0.31±0.01 


e 


1997.58 


765.0±38.3 


00.01±0.01 


Core 


1997.58 


164.0±8.2 


0. 1 1 ±0.01 40.7±2.0 0.09±0.01 


a 


1997.58 


22.0+1.1 


0.63±0.03 1.6+0.1 0.21±0.01 


b 


1997.58 


4±3 


1.60±0.08 11.6±0.60.01±0.01 


e 


1997.58 


4.0+0.8 


^ o 1 i 1 1 1 n 1 a 1 n A1 , n A1 

2.21±0.11 2.0±0.1 0.01±0.01 


xl 


2002.48 


1343+12 


00.01±0.01 


Core 


2002.48 


161+8 


0.16±0.01 35.6±0.20.09±0.01 


a 


2002.48 


18+4 


0.45±0.03 24.2±1.20.25±0.01 


b 


2002.48 


9.0+0.5 


1.23±0.04 19.0±0.8 0.29±0.01 


d 


2002.48 


3.0+0.5 


1.73±0.03 18.5+0.1 0.01±0.01 


e 


2002.48 


3.0+0.2 


2.52±0.03 22.2±4.9 0.24±0.01 


x2 


2003.20 


1342.0±6.7 


0.02±0.01 


Core 


2003.20 


159+8 


0.16±0.01 35.8±1.8 0.09±0.01 


a 


2003 20 


20+1 


0.45±0.02 24.3±1 .2 0.26±0.01 


b 


2003 20 


9.0+0.9 


1.23±0.06 18.3±0.9 0.31±0.02 


d 


2003 20 


3.0+0.3 


1.75±0.09 17.9±0.9 0.01±0.01 


e 


2003 20 


4.0+0.6 


2.76±0.14 9.7±0.5 0.16±0.01 


x3 


2003.49 


2819±2 


00.01±0.01 


Core 


2003.49 


128+1 


0.18±0.01 37.9±1.9 0.04±0.01 


a 


2003.49 


5.0+0.3 


0.47±0.02 21.4±1.0 0.01±0.01 


b 


2003.49 


8.0+0.4 


0.85±0.04 33.7±1.7 0.04±0.01 


c 


2003.49 


5.0+0.5 


1 .87±0.09 20.1 + 1 .0 0.25 ±0.01 


e 


2003.49 


5.0+0.8 


2.57±0.13 5.6±0.3 0.22±0.01 


x2 


2003 88 


2309+1 1 5 


0.04±0.01 


Core 


2003 88 


804+75 


0.17±0.01 30.5±1.6 0.10±0.01 


a 


2003 88 


14.0+0.5 


a 1 a a 1 ") 1 Oil /: a t 1 a r\o 

0.88±0.01 31.3+1.6 0.38±0.02 


c 


2004.60 


1444 ±72 


0.02±0.01 


Core 


2004 60 


328+16 


0.15±0.01 27.3±1.40.05±0.01 


a 


2004 60 


72+4 


0.53±0.03 24.6±1.20.21±0.01 


b 


2004 60 

_ v./ v./ r ■ v v./ 


30+2 


0.98±0.05 28.4±1.40.20±0.01 


c 


2004 60 


4.0+0.4 


1.8o±0.10 27. 7± 1.4 0.01 ±0.01 


e 


2004.60 


3.0+0.6 


2.86±0.1421.2±1.1 0.01±0.01 


x3 


2004.97 


907+45 


0.02±0.01 


Core 


2004.97 


475 ±24 


U.1U+U.U1 ZO. I±Y.D U.U0+U.U1 


a 


2004.97 


45 ±2 


M /1Q_i_M f\l 1 A -uC\ 1 M OO-lM M1 
1.4±U. 1 U.Zz±U.Ul 


D 


2004.97 


29±3 


1.08±0.05 25.1 ±1.3 0.35±0.02 


c 


2004.97 


4.0±0.4 


1.90±0.0921.5±1.10.01±0.01 


e 


2005.09 


376±20 


0.06±0.01 


core 


2005.09 


41.0±0.5 


0.26±0.01 8.5±0.20.16±0.01 


cl 


2005.09 


31+2 


0.91±0.03 20.5±0.20.64±0.03 


c 


43.2 2002.48 1484.0+74.2 


00.01±0.01 


Core 


2002.48 


329.0±16.5 0.040±0.00251.2±2.6 0.03±0.01 


a 


2002.48 


130.0±6.5 


0.18±0.0134.7±1.7 0.09±0.01 


b 


2002.48 


9.0±4.5 


0.61 ±0.03 22. 1 ±1 . 1 0.22±0.01 


d 


2002.48 


4.0±0.4 


1.14±0.06 7.1±0.40.10±0.01 


e 


2003.20 1853.0+92.7 


0.01±0.01 


Core 


2003.20 


211.0±10.60.080±0.00470.2±3.5 0.03±0.01 


a 


2003.20 


127.0±6.4 


0.15±0.01 38.4±1.90.13±0.01 


b 


2003.20 


6.0±0.3 


0.54±0.03 27.4±140.24±0.01 


d 


2003.20 


4.0±0.4 


1.08±0.05 31.1±1.60.01±0.01 


e 


2003.49 3310.0+16.6 


00.01±0.01 


Core 


2003.49 


294.0±14.70.080±0.00441.4±2.10.01±0.01 


a 


continued in next column 



S. Britzen: Non-radial motion in the TeV blazar S5 0716+714 



v Epoch 


c 

o 


r 9 M.A. Comp. 


[GHz] 


[mjy] 


[mas] [deg] [mas] 




2003.49 


86.0+4.3 


0.18+0.01 47.8+2.4 0.10+0.01 


b 


2003.49 


9.00+0.45 


a ii i a m a c a i i i a A1 .aai 

0.33+0.0245.9+2.3 0.01+0.01 


c 


2003.49 


4.0+0.2 


0.64+0.03 52.5+2.6 0.01+0.1 


d 


2003.88 


2072.0+10.4 


0.02±0.01 


Core 


2003.88 


1636.0+81.8 


0.040±0.00212.2±0.6 0.01 ±0.01 


a 


2003.88 


570.0+28.5 


0.13±0.07 29.4±1.5 0.09±0.01 


b 


2003.88 


57.0+2.85 


0.39±0.0229.9±1.5 0.26±0.01 


c 


2003.88 


4.0+0.2 


0.95±0.05 28.2±1 .4 0. 19±0.01 


e 


2004.601871.0+93.6 


0.02±0.01 


Core 


2004.60 


200±10 0.060±0.003 46.8±2.3 0.01 ±0.01 


a 


2004.60 


217.0+10.9 


0.16±0.0126.0±1.3 0.06±0.01 


b 


2004.60 


33.0+1.7 


0.36±0.02 12.3±6.20.20±0.01 


c 


2004.60 


22.0+1.1 


0.55±0.03 27.3±1.40.09±0.01 


d 


2004.60 


27.0+1.4 


0.83±0.0429.9±1.5 0.19±0.01 


X 


2004.60 


4.0+0.4 


1.07±0.05 25.5±1.3 0.01±0.01 


e 


2004.97 1382.0+69.1 


00.01±0.01 


Core 


2004.97 


197.0±9.9 


0.11±0.0123.4±1.20.08±0.01 


a 


2004.97 


12.0+0.6 


0.57±0.03 26.2±1.3 0. 10±0.01 


d 


2004.97 


18.0+1.8 


1.13±0.06 28.0±1.4 0.16±0.01 


e 


2005.09 


325.0±16.3 


0.05±0.01 


Core 


2005.09 


12.0+0.6 


0.33±0.0229.7±1.5 0.01±0.01 


c 



